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Abstract 

A simple model of a passive scalar quantity advected by a Gaussian non-solenoidal ("compress- 
ible") velocity field is considered. Large-order asymptotes of quantum-field expansions are inves- 
tigated by instanton approach. The existence of finite convergence radius of the series is proved, 
a position and a type of the corresponding singularity of the series in a regularization parameter 
e are determined. Anomalous exponents of the main contributions to the structural functions 
are resumed using new information about the series convergence and two known orders of the e 
expansion. 
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1 



I. INTRODUCTION 



The renormalization group can be considered as a most powerful method of investigation 
of the critical and scaling behavior. It produces results in a form of some, usually asymp- 
totic, expansions and only a few first terms are known analytically. Different resummation 
techniques are used to obtain reliable results [1]. The large-order asymptotic information 
of perturbation series for field-theoretic models is the base of critical exponents and scal- 
ing functions series resummation [1]. This information has been obtained with the aid of 
instanton approach at the example of the static 4 model and applied to the resummation 
problem [2, 3, 4]. 

Unlike the static models with 5-6 terms of a regular expansion for critical indexes known 
the perturbation series in the dynamical ones are calculated up to the third order only 
[5, 6, 7]. Nevertheless the resummation problem arises here also and demands an information 
about the large-order asymptotic in dynamic field theories. The aim of this paper is to 
develop the instanton approach for large-order asymptotic investigation in dynamics at the 
example of a simple dynamical model constructed from Langevin equation in the framework 
of the MSR formalism [8]. 

A model of the passive advection of a scalar admixture by a Gaussian short-correlated 
velocity field, introduced by Obukhov [9] and Kraichnan [10] attracts now a considerable in- 
terest. Some first structural and response functions in this model demonstrate an anomalous 
scaling behavior and the corresponding anomalous exponents can be calculated explicitly 
as within regular expansions in different small parameters as using numerical simulations. 
Thus this model provides a good testing ground for various concepts and methods of the 
turbulence theory like closure approximations, refined similarity relations, Monte-Carlo sim- 
ulations and renormalization group investigation. 

The advection of a passive scalar field </?(x, t) is described by a stochastic equation: 

d t (f - is Aip + gdi(viip) = f(x,t), (1.1) 

where dt = d/dt, di = d/chti, uq is a molecular diffusivity coefficient, A is Laplace operator 
with respect to x, £ is an artificial Gaussian scalar noise with zero mean and given correlator 
v(x, t) is a velocity field, g is a coupling constant; the sum in repeating indices here and 
henceforth are implied. Besides in the Kraichnan model the field v(x, t) obeys a Gaussian 
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distribution with zero average and the following correlator [11] 



D«(x - x') . (v ( (x, t )v j( x', t)) = DM - / (^ t'tlr^' '""'"'' 

where P J -(k), Pll(k) are projectors with respect to vector k direction, d is a dimension of 
the space, e is the regular expansion parameter and an ultraviolet (UV) regularizator of 
the model just as m is an infrared one; a is a coefficient representing non-solenoidal modes 
contribution, D is an arbitrary amplitude regular in e. The tensor indices of correlator will 
be omitted and implied henceforth. 

Using MSR formalism one can formulate the Kraichnan model as a field theory with the 
following renormalized action [11] 

S R = ^f- + - to - gdtvtp) + vZ v t^ - ^f^, 

integration in field arguments is implied. The renormalization constant Z v is determined by 
the requirement that the 1-irreducible part of the function (<p'<p) expressed in renormalized 
variables is UV finite (i.e., is finite at e — > 0). This requirement remains for Z u an UV finite 
arbitrariness; the latter is fixed by the choice of a renormalization scheme. The minimal 
substractions (MS) scheme show it's worth as the most useful for calculations. In this scheme 
all renormalization constants have the form "1 + only poles in e." The 1-irreducible part of 
the function (y?V) °f Kraichnan model is known exactly. It gives 



Z v = 1 — u _D C 



Sd d — 1 + a 

{1.6) 



=0 (2vr) d 2de ' 

where u = g 2 , S d = 2ir d / 2 /T(d/2) is the surface area of the unit sphere in <i-dimensional 
space. Note that the result (1.3) is exact, i.e., there are no corrections of the order u 2 ,u 3 , 
and so on [11]. 

It was stated recently [12, 13] the absence of an instanton within the MSR approach in 
Kraichnan model and it is the Lagrangian variables that can be used at the steepest-descent 
investigation of the structural functions < [<p(xi,ti) — (^(a^,^)]™ > in the large n limit. 
We deal with another problem here, namely the large-order asymptotic investigation for 
the perturbation series of the anomalous dimensions at arbitrary n. Nevertheless we also 
introduce the Lagrangian variables following [12, 13] to construct an instanton. 



Let us remind the main features of the Lagrangian variables approach. The Green func- 
tion G(x,s;y,t) corresponding to the Eq. (1.1) obviously satisfies 

(<9 S - vqA)G(x, s; y, t) + gdi(vi(x, s)G(x, s; y, t)) = 5(x - y)5(s - t). 
We can eliminate the (^-functions in the right-side by taking the following initial object 

G(x,s;y,t) = 0(s-t)P(x,s;y,*). (1-4) 

This yields 

{d s - z/ A x )P(x, s; y, t) + gd(y(x, s)P(x, s; y, *)) = 

P(x,s = *;y,*) = 5 d (x-y). 

The last equation may be considered as a Fokker-Plank equation of some different model 
with P(x, s) being a simultaneous distribution function. The term to term comparison 
allows us to reconstruct this model and to write down the solution P(x, s; y, t) explicitly: 

P(x, S ;y,*) = (5 d (x-X( S ))> c (1.5) 
d s X(s) = -gv(X(s),s) + ((s) (1.6) 
X| s= i = y, P c = 2is , 

< ... >^ denotes an average in the auxiliary random force ( with Gaussian distribution 
introduced is a corresponding correlator). Equations (1.5,1.6) describe particle motion 
in a random force field (. We see that the analysis can be reduced to the study of some 
interacting particles moving in a random medium as all the quantities in Kraichnan model 
can be expressed via Green function C(x, s;y, t) in the random field v. The medium is 
inhomogeneous: the Eq.(1.6) would be exactly the ordinary Brownian motion equation if 
there wasn't a term containing the velocity field that depends explicitly on the particle 
position. 

Let's insert in the Eq.(1.5) an additional path integration over an auxiliary field c(s) 
having in integrand Dirac 5 F - function with respect to the field variable c: 



P(x, s; y, t) = (J DcS F (c(s) - X(s))6 d (c(s) - x)> 
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(1.7) 



and rewrite 5 F - function in the form including a functional determinant: 



5 F (c(s) - X(s)) = det (d, + g^j 5 F (d s c + gv(c, s) - (). 

Then the d - dimensional 5 d - function in Eq. (1.7) can be transformed to the boundary 
conditions for c(s): 



P(x,s;y,t) = < J Dcdetld s + g-^\S F (d s c + gv(c,s)-C)) c 
c(t)=y ^ C ' 

Converting 5 F - function to a Fourier-like path integral form we get 



P(x,s;y,t) = M J DcDc' exp(-z/c' 2 + ic'c + igc'v(c, r)) det ( d s + g-^ J , (1.8) 
c(t)=y ^ C) 

the sums in the vector indexes of c(r), c'(r), v(c, r) fields and integration in field arguments 

are implied. Henceforth c denotes dc(r)/dr, M is a normalization factor appeared due 

to a functional determinant of the last 5 F - function transformation. The field theory 

obtained coincides exactly with the standard one. This can be proved by comparison of 

the diagrammatic expansions after the proper regularization of the functional determinants 

arisen. Namely det(<9 s + g5v/5c) can be considered as a constant independent on v, c fields, 

like in a usual MSR formalism [11] and it can be included in M factor. Thereupon M is 

determined by the free theory (at g = 0) and can be restored by comparison with it: 



(¥>(x,foV(y,t)) 



c(t )=x 



= M 

9=0 ciM y 



J DcDc' exp ^ - vc' 2 + ic'cj 



1 exp Y) ) , T = t -t. 



Usually the perturbation expansions series in the quantum-field theories turn out to be 
asymptotic ones with zero radius of convergence [1] . The situation is different in Kraichnan 
model: the exact result for ip 2 composite operator shows that these series have a finite radius 
of convergence [11]. In this work we prove the following iV-th term of the quantum- field 
expansions of Kraichnan model large N limit behavior 



(N) 



N b 



(1.9) 
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where b characterizes the type of singularity and a corresponds to the radius of convergence. 
Our aim is to calculate a and b coefficients for anomalous dimensions of composite operators 
(p n . The large-order information completes the perturbative one and allows us to extract 
singularities from perturbation series and to improve the convergence of the quantum-field 
expansions. 

The article [12] results in an interesting statement about the saturation of the scaling 
dimensions A n of the (f n operators. Namely they are predicted to tend to the finite value at 
n — > oo. One of our purposes is to verify this statement for the compressible velocity case 
on the basis of the results of the resummation for two orders of e expansion. 

The paper is organized as follows. In Section II the method of renormalization constants 
calculation is discussed. The general stationarity equations in field variables are obtained 
in Section III. The solution of these equations is described in Section IV. In Section V we 
deal with the stationarity equations in coupling constant and coordinate arguments. The 
integration over the scale parameter is described in Section VI. Section VII is devoted to a 
reexpansion of the anomalous dimensions of the composite operators ip n . The summary is 
written in Section " Conclusions" . 

II. LARGE-ORDER INVESTIGATION OF RENORMALIZATION CONSTANTS 

The quantity of interest for Kraichnan model is, in particular, the infrared behavior of 
the single-time structural functions 



It is determined by anomalous dimensions 7^ of the composite operators cp n , the former 
have been calculated up to two orders of e expansion [11]. 

Let us consider the response functions J dx dt ((p n (yi Q , to)v 9 '( x i^i) • • • f'i^n, t n )) repre- 
sented in a form 



<b(t,x)-v?(t,x')n, r=|x-x' 




(2.1) 



6 



Large-order investigation of this expression as a series in g is very difficult because of the 
presence in the Eq.(2.1) additional parameters such as m, x i? ti, e that can constitute different 
combinations comparable with the large parameter N of the steepest-descent method. The 
same problem for the classical tp 4 static model was described in [4]. Fortunately, as we are 
interested in the anomalous dimensions 7^n we can limit ourself by the consideration of 
renormalization constants Z v « of the composite operators ip n : 



•••</?')* = V<V?V--V>- (2.2) 

Constants Z^n are m and v independent in the MS scheme and contain only poles with 
respect to e variable. They are connected with the UV divergences of the equal time (ti = 
... = t n ) diagrams at zero external momenta. In the framework of e - regularization that 
is similar to the dimensional regularization so we can explore an analytical continuation to 
the e < region. In this region it is possible to set m = in the Eq.(1.2). This yields 



D v (x) = D 



2" £ r(-e/2) 



(4ir) d/2 r{d/2 + e/2){d + e) 



(d — 1 + a + e)5ij + e(a — 1) 



Ix| £ X 



e < 0. 



(2.3) 



An arbitrariness of the D parameter corresponds to the finite renormalization that does 
not affect the critical exponents. Let's then choose this parameter equal to 



_ T{d/2 + e/2)(d + e)¥ 
°~ r(l-e/2)r(d/2)d ' 1 ] 

so that the first fraction in the Eq.(2.3) contains a poles in e only. It is convenient to divide 
then the correlator discussed into three parts with different e - dependence: 



A(x) = 



D(x) 
2(d- I 



D v (x) = D(x) + A(x) + A(x), 

2|x| £ / , s x,Xi 



d(4n) d/2 T(d/2) 
a) 1 - |x| £ 



[a 



1 



d{47i) d/2 T{d/2) e 



A = 



2(d- 1 



"^_2 
x 2 

- a) 



d{An) d/2 T{d/2)e' 



(2.5) 
(2.6) 

(2.7) 
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Article [4] was devoted to the modification of the steepest descent approach for the 
large-order asymptotes investigation of renormalization constants at the example of the 
well-known static (/? 4 model in MS scheme and dimensional regularization. According to this 
article UV divergences have a sense as the perturbation objects and must be treated as the 
fore-exponential factors. So the factor exp(Sdiv) where Sd iv absorbs all divergent terms of 
the action has to be expanded in McLoran series in e before the application of the steepest 
descent approach. As the expressions (2.7) contain the UV divergences (at e — > 0, x — > 0) 
one need to include the corresponding terms to the divergent part of the action S^w and 
consider them as fore-exponential factors. In Section VI we will see that these terms affect 
only an amplitude of the asymptotic investigated. 

For f> 2 operator we can write in the Lagrangian variables with help of Exprs.(2. 1,1.4,1. 8) 



(V? 2 (x , t o y fa, ty (x 2 , f)) = M 2 J Dv J D Cl Dc 2 x 

Cl(t)=Xi,C2(t)=X2 

xDc[Dc 2 exp y - ^vD~\ - vZ u c' 2 - vZ v c 2 + ic[ci + ic' 2 c 2 + 

J Dvexpt-iv/^v) 



+igc' 1 v(c 1 ) + igc' 2 v(c 2 ) 



(2.8) 



c(t )=0 



• c(t)=0 

Integration of the Expr.(2.8) in v field yields 



(2.9) 



Cl(io)=C2(*o)=XO , 

M 2 J Dc 1 Dc 2 Dc' 1 Dc' 2 exp ( - vZ v (cf + c' 2 ) + 

Cl(t)=Xl,C 2 (t)=X 2 



+icj_ci + ic 2 c 2 — uc 1 



D(c 1 -c 2 ) + A(c 1 -c 2 ) + A 



c 2 . 



(2.10) 



It's worth while to scale u by v to obtain a new fully dimensionless coupling constant u. 
After the expanding of the factor Sdiv with UV divergent part of action 



Sdiv = -wciA(ci - c 2 )c 2 - uvc^Afci - c 2 )c 2 - v(Z v - l)(cf + c 2 ^ 



the Expr.(2.10) has a form 
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where 



Cl(to)=C 2 (to)=XO 

M 2 J Dc^Dc'^A^B^e-^, (2.11) 

Ci(t)=Xi,C2 (t)=X2 



= z/(cf + cf) - icici - ic' 2 c 2 + «i/ciD(ci - c 2 )c 2 , (2.12) 
A [ 2 ] = g (-wc^d - c 2 )4)J (2 13) 

3=0 ft 

Bm = g (-w.Ac,-^-!)^; 2 + c 2 2 ))^ (2 M) 

The similar expression can be written for ((/? n (x , t )(p' '(xi, i) . . . (p'(x. n , t)) Green function. 
It can be represented via the path integral in Cj, fields (i = 1, n) with the action 



sm = „ £ c f - i f; c& + ^ 2 < D ( C * - ( 2 - 15 ) 

i=l i=l iy^l 

In the last term the vectorial indexes of the fields and the tensor D are implied to be 
contracted. The normalization factor equals to M n now. The boundary conditions are 
c i(to) — x 0) c i{t) — x «- The pre-exponential factors and B'"' are similar to (2.13,2.14) 
with additional sums over % index of the fields Cj,c^. Note that Cj,c^ fields play a role of 
generalized coordinates and momenta for the action S^. So we can represent our system 
with a set of n moving quasi-particles. 

Now let's discuss the connection between the Green functions considered and the critical 
exponents. The singularities of the Green function and renormalization constants of the 
latter can be related via the identity (2.2). The renormalization constant Z^ determines 
uniquely the anomalous dimension 7^. 

The l.h.s. of the Eq.(2.2) is free of singularities and so does the right one: the singularities 
of Z^n are canceled by those of non-renormalized function ((p n (p' ...<//). We could use this 
to compute the singularities of Z^) by stating: 

• • • V') (0) + W • • • V?') (1) + • • • + 4°^V • • • y?'> (A °) = finite. (2.16) 
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Henceforth denotes the N-th order of the expansion for the value X in u. Usually 

expansions in the quantum-field theories are asymptotic ones with an exponential growth of 
the coefficients. Such growth would make all terms in the Eq.(2.16) except the first and the 
last ones irrelevant at N — > oo. Then Eq.(2.16) could be easily resolved for coefficients. 
The existence of the finite radius of convergence will be shown for Kraichnan model. Then a 
large set of terms in (2.16) become relevant and we face the problem of the unknown terms 
(ip n (p' . . .(p')^\i < N: we cannot pick out all the singularities present in the Eq. (2.16). 
That's why we shall use the identity 

ln(</?V • • • V') R = - In V + ln (v?V • • V) 

The l.h.s of the expression is finite again at e — > 0. To restore the N-th coefficient of the 
expansion for In Z^n in u we will investigate the residue in e for McLoran expansion of 
\n((p n (p' . . . (f'}. The "replica trick" that is based on the identity 

ln(y?V--V> = lim|-^V--V> L 

will be used to treat \n((p n (p' . . .ip'). The Green function in L - power in the last expression 
is then substituted by the path integral representation of (ip n ip' • • • V 9 ') with all integration 
variables considered now as L - dimensional ones. 

III. INSTANTON ANALYSIS AND IMMOVABLE PARTICLES. 

As in classical works [1, 2] we add an integration § du/u N+1 to produce the N-th order 
of the perturbation expansion. Summarizing all the remarks mentioned above the N-th 
order term of the expansion of In Z v n in u can be written as 

(In W ■ • V» ( ~> = ^ £n §~ L residue / ^ fj / «Di / <*x„„ x 

X J rf(x 2a - X lQ ) J rf(x 3Q - Xi a ) ...J d(x nQ - X lQ )M^ X 
t 0a )=x 0a } 

J [n^^< Q jANBNexp(-^] - ik(x 2a -x lQ )) (3.1) 



{Cj a (toa)=XOa} 

X 

{ c i.a (^q )— x iq: } 



with the action (2.15). The variables T, Xj, Cj, c- have now an additional replica index a = 
1...L. 
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Let us shift the variables to eliminate all dependencies of the path integral limits in 
coordinates 

/ \ - / \ _|_ ( x ia x 0o)(' r a t<x) /q <-,\ 

^ a 

the new fields c ict have zero boundary conditions. Sometimes we will return to the original 
notation c* and we will omit the index a for brevity. 

The steepest descent approach must be applied to the expression (3.1) with respect to 
all integrations except the one over the scale parameter of the model that has an essential 
non-instanton form. Indeed, overall UV divergences of interest arise from the integration 
over the scale parameter and can be extracted with the help of the integration by parts. We 
will show that at the instanton solution the value y = |x 2a — x lQ | doesn't depend actually 
on a and turn out to be the scale parameter. 

We proceed now to the finding of the saddle-point of the action in (3.1). Let us scale 
the variables to figure out the N dependence of the action. The renormalization constants 
considered are independent on renormalized diflusivity v and k momentum so they can be 
also scaled: 



c[ -> JVcJ k -> iVk, v = rj/N (3.3) 

(we perform this scaling both for the Expr. (3.1) and for the factor M integral representation 
(2.9); the Jacobians arising are constant and cancel mutually). Such a scaling produces 

N ndL/2 factor for the 

expression considered. The value has now the following form: 



S [ : ] = N K £ cf - i £ c'c, + y^2 <D(ci - c,)c; ee NSW , 
[ i=i i=i i^j J 

all the fields are assumed to have the replica index a. In the framework of the instanton 
approach we vary the functional 



S = N E (s* 1 + M*2a ~ x la )) + N In 



u 



with respect to all the variables. The variations in c m , c' rn yield 



11 



— -0 =► it' -gnVc' dD(c m - Cl ) 

Sc m -° * ^m-uv^c m d{Cm _ Ci) c„ (3.4) 

= =>• ic m = u?7^^(c m - c/)cj + 2?7c^. (3.5) 

The dependence on coordinates Xj is represented in the functional S only by the fields 
Cj, Cj (see the Expr. (3.2)). So the variation in x yields the full momentum conservation 
law for the quasi-particles: 



0x , =1 
The variations with respect to the coordinates x m — xi (m > 3) yield: 



. to to 



- J drc m + u V J dr^r- ^ c m _ £ = 0, m > 3. (3.6) 

1 t 

Combining together Eq.(3.4) and Eq.(3.6) we get the equation 

1 t to 

f J drc' m --/ rfr(t - r)tf m = <,(*) = 0. (3.7) 
t t 

The last expression can be considered as a boundary condition at r = t on the field c^(r) 
(m > 3). The Eq.(3.4) is a first order differential equation with respect to c' m and has a zero 
boundary condition (3.7). It has then the trivial solution 

c' m ( T )=0, t<r<t m > 3, a = l...L (3.8) 
that is locally unique. The Eq.(3.5) provides: 



c m = c m (r) = x m = x m > 3, a = 1 . . . L. (3.9) 

This means that the instanton sought for implies only a couple of quasi-particles in 
motion while the others rest at the point x . The saddle-point equations reduce then to the 
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n = 2 case with only four nontrivial fields Ci,*^,^,^ (replica index a is assumed). Such 
a skewness of the instanton is explained by the term k(x 2a — Xi a ) of the variable action. 
Indeed, as we are interested in renormalization constants that don't depend on a conjugate 
momentum running through diagrams of the response function considered we have set all 
momenta except k in the Expr. (3.1) equal to zero. The momentum k can't be equal to 
zero since the object considered would turn out to be divergent. Note that this divergence 
is unrelated to the problem as it has a trivial power form and doesn't affect on the radius of 
convergence analyzed. Nevertheless the accurate treatment requires k ^ that causes the 
skewness obtained. 

Thus the analysis of the stationarity equations demonstrates the trivial solution (3.8,3.9) 
for the variables x^, k > 3. The problem reduces then to the y? 2 case and we can limit 
ourselves to the study of the saddle-point of the following nontrivial integral representation 

_1^_^_]7_ J rf Xo | dxJdT J Dc 1 Dc 2 Dc , 1 Dc' 2 A^B^M 2 e- s{ "\ (3.10) 
m u a ' o 

that represents the case n = 2. Indeed, the substitution of (3.4,3.5) to the expressions for 

reduces their values to the two-particle case A^ 2 \B^ and S$ as well. 

Let us introduce new variables: p = c[ — c' 2 , P = + c' 2 , q = Ci — c 2 , Q = Ci + c 2 . 

The fields q, Q satisfy zero boundary conditions. Sometimes we will use the notation q = 

c i _ c 2, Q = c i + c 2 for brevity. The variables introduced correspond to the center mass 

frame of reference of the quasi-particles. The expression (3.10) transforms into: 

J_ j J dXo y d^JdTM 2 A^B^ J DpDqDPDQe- N ", (3.11) 

o 

L 

S=Y,( S *+ ikx *)+\nu, x = x 2 -xi, (3.12) 



a=l 



s? = -ip(4 + g - ^p(q + 2X0 "^" X2 ) + u i (P + pMq)(p - p) (3.i3) 

All the variables in the r.h.s. of the Expr. (3. 13) are assumed to have the replica index a. 
Zero conditions on variations of the action S in q, Q, p, P, x, x and u give the instanton 
equations. 

It is convenient to carry out the calculation in two steps. First we solve the instanton 
equations in fields P, Q, p, q and the variable x . We use their solution to simplify the 
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functional S. Then we vary the simplified functional S in x a and u we obtain saddle-point 
equations in a simplified manner and we turn out to be successful in solving them. 



IV. TWO QUASI-PARTICLE SADDLE-POINT SOLUTION 

The saddle-point of the action S with respect to x , Q, P, q, p is determined by the 
following equations: 



SS_ 

5S 
6Q 
5S_ 
5P 
6S_ 
5q 



= 



IT- = =>• 



5p 



uri d 



p = o, 
p = o, 

iQ = wr/D(q)P + 2r/P, 
pD(q)p - PD(q)P 

iq = 2r/p — wr/D(q)p 



with the following boundary conditions 



q(t) = -x, Q(t) = xi + x 2 , 
q(*o) = 0, Q(t ) = 2x . 

To find the solution we suppose that the vector p is parallel to q. In fact one can show 
that this is the only solution of the original system satisfying the given boundary conditions. 
Similar to [12] the simplified system is trivially integrated: 



Q(r) = 2x = x 1 +x 2 , P(r) = 0, (4.1) 



P(r) = 2 ,4g q(r)) . qM = ^A-^(^)), (4-2) 

J^x) = f . dZ D(z) = ^ , (4.3) 

A. Pn-uriD{z) d(4rc) d/2 T(d/2y 

where D(z) and x is the projections of the tensor D(q) and of the vector x to the direction 
of the vector p. As usual the index a is implied for the time variables T, r, all fields and 
the variables x, x. 
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The functional S at the saddle-point equals 



^E&M-itoO+lnu, SM = 3^, (4.4) 

a=l 

where the replica index is indicated explicitly. 

Let's note that due to P = the following identity holds = — c' 2 and the A - term in 
B^ factor (2.11) is canceled by the renormalization (Z v — 1) of the diffusivity v and thus 
= 1 (2.14,2.7,1.3). The pre-exponential factor A® in the saddle point is given by 



g^ ^W/.M V J- A(z)dz 

V. A SADDLE-POINT SOLUTION OVER u AND x Q 

It is obvious that the variations in u and x Q of S given by the Expr. (3.12) and S given 
by the Expr. (4.4) are explicitly the same. The form (4.4) is much more convenient: the 
variation with respect to u, x Q 



— = o y ^ / i( x °) / 2(x a ) = l , 51 , 



dzD(z) 



<* x * 2T a ^2r] - uD(* a ) 



ik. (5.2) 



These equations include T a parameter: let's suppose their solution exist. Substituting it 
into the expression (4.4) we obtain the action S$ = NIf/T a i.e. the integration over T a has 
the form 

r°° ( NI? \ 

J Q dT a f(T a )e W l-^J. 

The convergence of this integral is assured by f(T) factor that absorbs all the fluctuation 
integrals over the fields and w,x ,x variables. It is easy to see that the main contribution 
to the integration over T a at iV — > oo is determined by large T a region. Let's return to the 
Eq.(5.1),(5.2) analysis. The l.h.s. of the Eq.(5.2) rests non-zero in the T a — > oo limit only if 
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the value y2r] — uD(x. a ) is small (~ T~ v ). On the other hand the Eq.(5.1) is fulfilled only 
for large values of I 2 , i.e as the denominator ^2<q — uD(z) reduces to zero at the point that 
tends to x + at large T a . Therefore the solution of the system (5.1), (5.2) can be written 
at large T a as 



S a 2 
X « =V + Tr U = D(y + SoEaV) (5 ' 3) 

with 5o, S a being some constants and y being an arbitrary scale parameter. The substitution 
of this solution into the Eqs.(5.1), (5.2) allows to calculate S ,5 a constants. In the limit 
T — > oo this saddle-point expressed in x a ,tt is 

k o 

x a = y-rz-r, a = 1 . . .L, u 



|k|' ••• ' D( y y 

The arbitrary scale parameter appearance is the usual fact at instanton investigation of scale 
invariant models [2]. Its correct consideration demands an insertion of a unit decomposition 
in the initial expression (3.1) 



/*oo 

1 = / dy5(y - |x Q= i|). 
Jo 



The integration over the scale parameter can be made explicit by Faddeev-Popov method. 
It turns out to have the non-saddle-point form and its singularities at e — > determine the 
renormalization constants and the critical exponents (the same situation is observed for ip A 
model [4]). The 5 - function mentioned above plays an important role as it solves the so called 
"zero mode" problem [1] in the carrying out of the integration over the fluctuations. Besides 
it produces an additional y/N factor as a result of the elimination of the integration over one 
of the degrees of freedom ("zero mode") in the numerator of the fluctuation integral. The 
other factor ]Sf-( ndL + 1 )/ 2 appears from the fluctuation integration over the variables u, x , x ia 
(i = 2 ... n, a = 1 . . . L). Together with N ndL l 2 factor appeared from the scaling (3.3) the 
power function in iV reduces to unit. 

The last question we have to discuss in this Section is the saddle-point method applica- 
bility. It should be mentioned that the fluctuation integration can be carried out at e = 0. 
Though the longitudinal projector 

(Cj Cj) a (Cj ) {, 



p, 
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in the correlator (1.2,2.3) cannot be expanded in the series of the fluctuations of immovable 
quasi-particles variables (i, j > 2) with respect to the saddle-point, the integration over 
these fluctuations (Sci, Sc'^i > 2) is Gaussian and can be correctly performed by the means 
of transition from the path integral in <5c^ (i > 2) to the integration over 5c^i,5c^ ± fields 
that represent the longitudinal and transversal projections of <5c^ field on the Sc^ direction. 

The consistent analysis of the higher variations of the action (2.15) around the saddle 
point demands their examination in the large T limit as it is this region that contributes to 
the integration over T of the expression studied. Using the explicit solution (4.1,4.2) one 
can show that the higher variations of the action considered rest finite in the large T a limit. 
This assures the consistency of the instanton approach and the existence of the fluctuation 
integral. 

VI. THE INTEGRATION OVER THE SCALE PARAMETER 

Finally the result of instanton approach is 



The factor K arises from the L - dimensional integration over T. The function 
n(y)(n(0) < oo) results from all the fluctuation integrations, besides the factor Ila 6 *' 1 ^ is 
also included in k. The factor 1/y is extracted from fluctuation integrations without explicit 
calculation by the means of the dimension analysis. Indeed, the value lnZ^n is determined 
by logarithmic divergences of the diagrams corresponding to (ip 2 ip' • • • f')^ function. This 
gives the logarithmic behavior of the integral in y (6.1). The convergence of the Expr. (6.1) 
at large y is ensured by re(y) factor. Let's remind that this factor also depends on L. 

Using the formula (5.3) we note that /a(x q ) ~ T a at the saddle-point, so the Expr. (4.5) 
can be rewritten as 



with K(e) and f(e) being regular functions in e, K(0) ^ 0, /(0) = 1. The extraction of 
UV divergences (i.e. poles in e) have to be executed in the framework of the formalism 
developed in [4]. The parameter Ne is supposed to be small. Using the explicit formula for 
D(y) (4.2) the divergences at e — > can be easily extracted by integration by parts: 



In 



2$ ~ hm A residue \ K L f d M^)\(y) f[ 



i 



A P] 



(6.1) 
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o y V {y) ei ~e^h(N + s)s\{j-s)V ^ 
Due to the choice of D v correlator (see definition of D constant by the Expr. (2.4)) 
the studied object (6.1) doesn't contain Ne corrections. Then term doesn't affect the 
position and the type of the singularity and contributes at large order in 1/N only to the 
amplitude of the asymptotes. The final result for the single pole in e of the renormalization 
constant series in u is given then by the following formula 



a 

(An) a/ *T(d/2)d , 



N 



VII. ANOMALOUS DIMENSIONS RESUMMATION 

We proceed now to the critical exponents A^n analysis. The general expression for critical 
exponents is 



A^n = dpn + = 7l(-l + S) + 7^" , 

where d^n and 7^ stand for the canonical and anomalous dimensions of the composite 
operator y? n . Our main goal is to analyze the expansion of the anomalous dimension 7^ 
in e using the information on the series singularities. These singularities are of the main 
interest since they can be used for resummation of the series. General theory [1] provides 



d 

7 v n = D^hiZ^n = (3{u)— residue [InZ, 



d 

—2u — residue [In 

du e=Q 



with u* being a fixed point. The latter is known exactly for Kraichnan model: 



(7.1) 



The substitution of the Expr. (6.3, 7.2) to the formula (7.1) gives the transformation of u 
- expansion of the anomalous dimension 7^ to the series in e considered at the large order 
of perturbation theory: 
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1^ 



(TV) 



JV>0 

-2a 



-i N 



(d-l + a) 

The absence of a power term iV b in the resulting asymptote (i.e. b = in the Expr. 
(1.9)) indicates the simple pole singularity. The convergence radius R c and the position of 
the nearest singularity e c for the anomalous dimension expansion in e is easily found to be: 

d — 1 + a 



2a 



This information is used for resummation the 7^ series: we extract the singularity of the 
expansion adopting simple rational representation: 



1^ 



S ~(fc) k 

<p* £ - e e ■ 



(7.3) 



fc=i 



Term by term comparison of these expansions for 7^ yields the following identities for the 
few first coefficients 7^ and 7^ 



ly> n — /^"fcc; lip™ — l^™ 7^ nfc c; 

We now need the values of the coefficients 7^ to carry out the resummation. The anoma- 
lous dimension 7 v n,n > 1 is calculated up to the second order in [11] and the coefficients 
are given by the formulae: 



(1) an(n — l)d 
a — 1 + a 

(2) a(a — l)n(n — l)(d — 1) a 2 n(n — l)(n — 2)dh(d) 
lipn = (d-l + a) 2 + (d-l + a) 2 ' 

w = g iH1+im :.. (k+m =w+ «</* m, 

where F represents generalized hyper-geometric function. The resummation is then straight- 
forward and gives the following coefficients 7^, k — 1,2: 
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7'n = -n(n - 1)-, 



lip™ 



— n(n — 1) + 



a(n — 2)dh(d) 
2(d-l + a) ' 



and the anomalous dimensions equals to: 



= 



2an(n — \)e 
' d- l + a(l + 2e) 



d 



+ e 



ad(n — 2)h(d) 



e + 0{e 2 ) 



(7.4) 



2(d- 1 +a) 

Fortunately, this result can be verified: for n = 2 an exact value of the anomalous 
dimension is known: 



7^2 = 



2ea(d + 2e) 



d- l + a(l + 2e)' 
while the resummation (7.4) for anomalous dimension 7^2,72 = 2 yields 



7^2 



2ae [d + 2e + 0(e 2 



d-l + a(l + 2e) 

We see that no high order corrections to the resummation result occur in this case; resum- 
mation gives the exact answer. The singularity found for Z v 2 turned out to be unique. This 
case illustrates the information about the type and the position of the singularity to allow 
us appreciably supplement the information obtained from direct perturbation expansion. 



VIII. CONCLUSIONS 



This article demonstrates that the instanton approach applies to the large-order analysis 
of the dynamical model. We can state that the study of the behavior of the convergent series 
is a more difficult problem than that for the divergent ones. Using the instanton approach one 
can state that the perturbation series for the scaling dimensions of the operators considered 
are convergent with the finite convergence radius calculated. In each case taking into account 
the information about singularity type of the series permit to accelerate considerably the 
convergence of these series. 
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Concerning the saturation problem [12] one can say that our results don't prove its 
existence. It is possible that two orders of the e expansion are not sufficient to observe the 
saturation. Another possibility is that there is no saturation in the studied model (saturation 
was predicted for the Kraichnan model with the transversal velocity field a = [12], where 
our results can are not valid). 
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